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Abstract 

Inspired by the warped Randall Sundrum scenario proposed to solve the mass scale 
hierarchy problem with a compactified fifth extra dimension, a similar model with 
no metric singularities has been elaborated. In this framework, the Kaluza-Klein 
reduction equations for a real massless scalar field propagating in the bulk have 
been studied carefully from the point of view of hermiticity so as to formulate 
in a mathematically rigorous way all the possible boundary conditions and cor- 
responding mass eigenvalue towers and tachyon states. The physical masses as 
observable in our four-dimensional brane are deduced from these mass eigenvalues 
depending on the location of the brane on the extra dimension axis. Examples of 
mass towers and tachyons and related field probability densities are presented from 
numerical computations performed for some arbitrary choices of the parameters of 
the model. 
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1 Introduction 



In a previous paper [T], we reanalysed mathematically, within the Arkhani- 
Ahmed, Dimopoulos, Dvali [2j large extra dimensions model, the procedure of 
generation of the Kaluza-Klein masses [3] , stressing that it is the momentum 
squared in the extra dimensions (and not the momentum itself) which is the 
physically relevant quantity and hence corresponds to an operator which must 
essentially be hermitian. For illustration purpose, we restricted ourselves to 
the case of a five- dimensional massless real scalar field supposed to propagate 
in the flat five-dimensional bulk. The extra dimension is compactified to a 
finite range, say [0, 2nR], either on a circle (then R is interpreted as the radius 
of the circle) or on a finite strip (then 2ttR = L is the length of the strip). 
All the allowed boundary conditions resulting from the requirement that the 
extra dimension momentum squared must be a mathematically precisely- 
defined symmetric operator have been established (see also [1]). We deduced 
from them, besides the usual regularly spaced Kaluza-Klein mass towers, new 
towers with non regular mass spacing and tachyons. These considerations 
should be extended to vector and tensor fields. 

In this article, inspired by the Randall and Sundrum scenario we have 
developed a model based on a warped space with one extra dimension and 
basic parameters chosen so as to solve similarly the mass scale hierarchy 
problem. 

In Randall-Sundrum, the fifth dimension s is compactified to an orbifold 
of radius R. A so called Planck brane is located at s = while the TeV brane 
or Standard Model brane is at s = ttR. We depart from the original scenario, 
postulating that the compactification is on a strip, that the metric has no 
singularity and that only one particular brane, the TeV brane on which we 
live as a four dimensional observer, has to be considered physically. 

In this framework and again restricting to a real massless scalar field 
propagating in the bulk we have carefully studied the hermiticity properties 
of the operators in the Kaluza-Klein reduction equations for the adopted five 
dimensional warped space. We have enumerated all the allowed boundary 
conditions and from them we have deduced the corresponding Kaluza-Klein 
mass eigenstate towers and tachyon states and have studied their main prop- 
erties. 

As will be shown, the values of the observable physical masses in the 
Kaluza-Klein towers can be deduced from the mass eigenvalues and depend 
on the particular location of our four dimensional brane in s as do the field 
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probability densities which are related (when dynamics and kinematics are 
included) to the overall probability that the associated mass states would 
appear to the observer. 

In the five dimensional bulk, we also postulate that all the dimensionfull 
parameters are scaled with a unique mass, the Planck mass Mpi. It then 
happens that, by an adequate choice of the reduced parameters defining the 
model, all the low lying physical masses obtained from the eigenvalues of the 
Kaluza-Klein reduction equation are of order TeV for an observer living in 
our 4-dimensional brane. 



2 Operators in the five- dimensional warped 
space. Mathematical considerations 

The warped five dimensional space with coordinates x A ,A = 0,1,2,3,5 is 
composed of a flat 50(1,3) invariant infinite four-dimensional subspace la- 
beled by x^ (/i = 0, 1, 2, 3) with signature (+, — , — , — ) and a spacelike fifth 
dimension with coordinate x 5 = s on the finite strip < s < 2ttR. The 
metric, unique up to rescaling, 

dS 2 = g AB dx A dx B = e~ 2ks dx^dx^ - ds 2 (1) 

satisfies Einstein's equations with a stress-energy tensor identically zero and 
a bulk negative cosmological constant A as the unique origin of the induced 
Riemann metric ([1]). Indeed k, here chosen positive (see the discussion about 
the sign of k in Appendix(|B])), and A are related by 



k = ^~, k>0. (2) 

A free massless scalar field <&(x, s) in this warped space satisfies the in- 
variant equation 

□Ricma„n$ = -^-d A ^gg AB d B ^ = . (3) 

From the metric ([1]), ^fg = e~ 4ks and Eq.Q becomes 

(^e 2ks U 4 - e 4ks d s e- 4ks d s ^j $(ac", s) = (4) 
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where = d^d^ is the usual four dimensional d'Alembertian operator . 

We now carry on a careful study of the hermiticity properties of the 
operators appearing in (j3j). For scalar fields, the invariant scalar product is 

/+oo /•2ttR 
<fx j ds y/gV*(x,s)$(x,s) . (5) 
-co JO 

Remember that an operator A with dense domain D(A) 

• is symmetric for a scalar product if 

= (6) 

for all the vectors G D(A) and <E> G -D(A), i.e. if the adjoint operator 
A^ of the operator A is an extension of A: A^<& = v4$ for all $ G -D(v4) 
and D(At) D D(A), 

• is self-adjoint if A^§ = A$ for all $ G -D(A) and moreover D(A^) = 
D(A), i.e. if the operator is symmetric and if the equation ([6]) cannot 
be extended naturally to vectors \I/ outside D(A). 

We will call a differential operator which is symmetric up to boundary 
conditions a formally symmetric operator. 



Symmetric Operators 

For the scalar product (jSJ), the operator □Ri ema nn © is formally symmetric. 
The two operators which appear in (j3j) have the following properties 



A x = e 2ks D 



Ao 



e 4ks d s e~ 



i 

4k S Q 



is self— adjoint 

is formally symmetric . 



(7) 



By partial integration, one finds that the domain condition for the formally 
symmetric operator A2 to be symmetric is 



HirkR 



2ttR 



(8) 



which means that any field $ and its derivative, both evaluated at s = and 
s = 2nR, have to satisfy the same specific linear relations. These relations 
which express the boundary conditions will be studied carefully later. 
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Commuting operators 

The two operators defined above (j7|) do not commute and hence cannot be 
diagonalized together. 

Multiplying equation on the left by e~ 2ks leads to the following equiv- 
alent equation 

- e 2ks d s e~ Aks d)j $(a/*, s) = (9) 

which defines two commuting operators with the following properties 

B\ = is self— adjoint 

B2 = e 2ks d s e~ iks d s is not even formally symmetric . (10) 



Puzzle 

So we are facing a puzzle 

• either the two operators Ai,A% ([7]) are formally symmetric but do not 
commute, 

• or the two operators Bi, B2 ffTOl) commute but the second is not formally 
symmetric. 

Solving the puzzle 

The puzzle can be solved remembering recent discussions about non-hermitian 
operators having real eigenvalues [6]. It was shown in [7J that in many cases 
these non hermitian operators are in fact equivalent to hermitian operators 
by a non-unitary change of basis. This is the case here. Indeed after some 
algebra, considering the non-unitary transformation induced by V(s) = e ks , 
one finds 

B x = VB X V- V = B 1 
§2 = VB2V- 1 

$ = 1/$ (11) 

with the result that Bi and B 2 are at least formally symmetric operators for 
the induced scalar product 

/+00 p2ttR 
d A x / ds e~ 6ks **(x,s)$(x,s) . (12) 
■00 Jo 
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It is also easy to see that the natural conditions for the operator B 2 to be 
symmetric when acting on the space of vectors $ are 

These conditions turn out to be fully compatible with those (jHJ) obtained for 
s) from the requirement that A 2 is symmetric. 
Thus even though the operator B 2 is not symmetric, it is equivalent 
through a change of basis to a symmetric operator B 2 and hence will produce 
real eigenvalues which will be related to the Kaluza-Klein tower masses as 
will appear below. 



A2nkR 



(13) 



3 The Kaluza-Klein reduction equations 

To simplify the discussion, we concentrate on the particular case of a real 
massless scalar field. The general procedure to solve the basic equation ([9]) 
along the Kaluza-Klein reduction method is well-known. One supposes that 
the field s) is a linear combination of terms where the variables x M and 

s separate 

<iK^s) = X>L*V)0i?OO • (i4) 

n 

Then s) is a solution of (JSJ if 

B x ^{x^)= D^V) =-m^£V) (15) 
B 2 $\s) = e 2fes 9 s e- 4fcs 9 s 0k s] ( S ) = -m 2 0W( a ) . (16) 

From the arguments given above one then concludes 

1. The operators B\ and B 2 commute and can indeed be diagonalized 
simultaneously 

2. The operator B 2 is equivalent to a formally symmetric operator through 
a non-unitary change of basis and hence, taking into account boundary 
conditions compatible with (El), Eq. (Tl6l) gets real m 2 eigenvalues which 
can be positive, zero or negative. 
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3. The operator Bi is self-adjoint. By (TT5|) . the solutions for > 
correspond to four- dimensional physical particles, those with m? n = 
to four-dimensional massless particles and those with w? n < to four- 
dimensional tachyons. 

The boundary restrictions (jHJ) can be conveniently rewritten 



When k = 0, these restrictions are identical to the restrictions applicable in 
the fully flat case which we studied in pQ. As a consequence, the boundary 
conditions compatible with (fT7|) can be obtained by simply replacing 4>n\o) 
and <9s0M(O) by respectively e 47rfcR 0n'(O) and e AnkR d s <p^(0) in the boundary 
conditions which we listed for the flat case. 

In fact, each field must satisfy the same set of boundary conditions which 
consist of at least two linear relations. This set defines specific domains in 
the Hilbert space. In Tabled!]), we give all the possible independent sets of 
boundary conditions expressed by just two linear relations. The box condi- 
tion appears in Case A6. 

Boundary conditions expressed with more than two linear relations can 
be considered as restrictions applied to the cases with two relations. The 
domain of the operator is then reduced. For three linear relations, the sets 
are given in Table([2]). The set involving four relations consists of <pn\o) = 
0, <p l n ] (2nR) = 0, d s <p l n ] (Q) = 0, d s <p l n ] (2nR) = 0, which can only be satisfied 
by the trivial field (pn\s) = and can thus be forgotten. 

The boundary conditions have to be imposed to the general solutions of 
the equation (fl6l) which are 

• For m\ > 0, the solutions are linear superpositions of the Bessel func- 
tions J2 and Y 2 (see Appendix (TAj) ) 



(#(^* 1 )-(^4 s1 )* 1 ) 



J 2nR 



= ((e A * kR ^)(d s e^ kR <P [ : ] ) - (d s e 47TkR ^ p s] )(e 47rkR <P [ : ] )^ .(17) 



J 
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• The solution for m„ = is 



^(s) = a e iks + r . (19) 

• The tachyon solutions for mf = —hf < are linear superpositions of 
the modified Bessel functions I 2 and K 2 (see Appendix (TAl) ) 

= (20) 

• In the above formulae, a n , r n , a , r , a t and r t are constants. 

4 Physical considerations. Masses 

We now extend our considerations to potentially physical consequences in- 
cluding for example the possible discovery of TeV warped states at high 
energy colliders. In the next subsections, we first discuss the magnitude of 
the parameters k and R which occur in the model and also the magnitude 
of the parameters ccj ... which define the boundary conditions, postulating 
that there is only one scale in the model, the Planck mass. We then write 
explicitly the equations determining the mass eigenvalues. We discuss the in- 
terpretation of these mass eigenvalues in terms of the physical masses as they 
would be observed in our brane. In particular the conditions for the existence 
of zero mass states and of tachyons mass states are deduced. Examples are 
finally given and discussed. 

4.1 The parameters 

The general philosophy underlying the warped approach, which was proposed 
to solve the hierarchy problem, is that there is only one mass scale, the Planck 
mass Mpi pa 1.22 10 16 TeV, and hence that any dimensionfull parameter p 
with energy dimension d is of order Mp X . More precisely 

p = p(M Pl ) d (21) 

with p a pure number of order one. In particular, k = kMpi and R = RMp^. 
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As stated above ([2]), the parameter k is chosen positive. As will be seen 
hereafter, it appears that for an observer sitting at s = a reasonable choice 
for the product kR is 

kR = kR w 6.3 (22) 

as the resulting lowest masses in the Kaluza-Klein towers would then be of 
order 1 TeV. 

4.2 The solutions 

We restrict ourselves to all the boundary conditions which are expressed by 
two relations (see Table (CQ)). We first discuss the Kaluza-Klein mass towers 
in general, then in particular the towers which have a mass zero as their 
lowest mass state and finally the towers with a tachyon state. 

Towers m 2 n > 

Two boundary relations being applied to the fields (fT8l) lead to two linear 
homogeneous equations in terms of the parameters a n and r n . The coefficients 
turn out to be linear combinations of Bessel functions evaluated for s = 
and s = 2ttR, i.e. with arguments respectively equal to 



F 2 = (23) 
In order to find non trivial solutions for a n and r n and hence non trivial 

[si 

fields (f) n , the relevant determinant has to be equal to zero. This leads to 
the mass equation whose solutions provide for chosen parameters fixing the 
boundary conditions the mass eigenvalues m n building up the related tower. 
Once these mass eigenvalues are known, the corresponding a n /r n ratios are 
deduced from one of the two original equations. For each set of allowed 
boundary conditions, both the mass equation and a cr n /r n relation are given 
in Table®. 

Zero mass states 

With two boundary relations being applied to the fields (fl9|) . one obtains 
again two linear homogeneous equations in terms of the parameters a and r . 
The condition that the relevant determinant is zero is in fact the constraint 
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that has to be satisfied by the boundary condition parameters for a zero mass 
state to exist. One of the two equations fixes again the ratio <Jo/ T o- 

For each of the allowed set of boundary conditions, both the parameter 
constraint and a <Jo/t equation are given in Table®. Since e 4nkR takes such 
a large value and since the reduced parameters di, . . . are assumed to be 
of order one, approximate relations, valid to a high degree of precision, are 
easily obtained. They are also listed in Table®. We should remark that 
there is no zero mass state for the box boundary condition (Case A6). 

The parameter constraint equation for a zero mass state defines a surface 
in the parameter space (see approximate equation in Table®). In many 
cases, if one follows a path in the parameter space which crosses the parame- 
ter constraint surface, the lowest mass eigenvalue goes smoothly toward zero, 
takes the value zero as the path goes through the surface and emerges as a 
tachyon state with low h 2 = — m 2 . Examples will be given in subsection f)4.4l) 
for the boundary Cases Al, A4 and A3. A different behavior shows up in the 
Case A3 when in particular pi passes zero. The zero mass appears suddenly 
as the surface is crossed. No small mass, no tachyons appear on either sides 
of the surface. This results from the fact that the solution m 2 = is of higher 
order. 

Tachyon states h 2 = —m 2 > 

The equations for h 2 can simply be obtained from the equations giving the 
tower masses by replacing the Bessel functions J n by I n and Y n by (— l) n K n . 
They are summarized in Table® for each set of boundary conditions. 



4.3 Physical interpretation of the mass eigenvalues 

At this stage of the discussion, in order to make connection with physics, one 
has to take into account the position s = s (0 < Sq < 2ttR) where the TeV 
brane, the brane in which we live, is located. Indeed the deduction of the 
physical masses in terms of the mass eigenvalues depends crucially on this 
position. 

In our brane, the space time part of the metric ® has a factor e~ 2ks °. 
By a change of variable 

x^ = e~ ks °x, (24) 
the space time metric in normal local units becomes dS 2 = dx^dx^. Hence 
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the physical mass mn, s obeys the equation 

□ 4 0W(F) = -ml tS0 JH(2?*) . (25) 
Comparing (I2"5"j) and (IT5"|) we find 

"Vso = efcs ° m n • (26) 

This gives the relation between the m n eigenvalues which appear in Eqs. (ll5j) 
and f)16p and the observable physical masses m njS0 in the brane. In the case 
s = 0, the eigenvalues and the physical masses are equal (m n> o = m n ). The 
low lying masses are of order 1 TeV when kR = 6.3. If so differs appreciably 
from zero one sees that m„ )S0 may become large enough to spoil the hierarchy 
solution. However, as will be seen later, the solution of the problem is fully 
restored by an adequate increase of kR 

4.4 Examples 

In this subsection, we show for illustration the results of some numerical 
computations for Kaluza-Klein mass eigenvalue towers and field probability 
densities corresponding to two sets of boundary conditions belonging to Case 
A4 and Case Al respectively. We also make at the end some introductory 
comments about the mass tower structure in the Case A3 

In all the examples, the parameter kR has been set equal to value 6.3 and 
the parameter k has been chosen equal to one for convenience, so 

kR = 6.3 , k = 1 . (27) 

The extension to other values of kR and of k is outlined. 

As stated in fj4.3j) . the computed mass eigenvalues would be the physical 
masses for a four- dimensional observer at So = 0. The physical masses for 
s 7^ can be deduced from (T2"oT) . 

Case A4 mass eigenvalues 

The case A4 is simpler since there is only one free parameter fixing the bound- 
ary condition. We remark that the mass equation for the tower eigenvalues 
(see Tablefll])) is invariant under the following rescaling with the arbitrary 
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parameter A 



R -> 

k -> 
k — ► 
m — > 

which leaves fci? = kR invariant. Hence 



R 

A_ 
Afc 

Ak 

Am (28) 



Xmi(k, k) = rrii(Xk, Xk) (29) 

which allows one to determine the tower mass states for other values of k 
from the mass states corresponding to our choice k = 1. Indeed 

mi(k,R) = krrii (l, =J . (30) 

The ten first mass eigenvalues mj (i — 1, . . . , 10) in the towers are given in 
the Table© (remember fT2~T|) ) and for a few values of the boundary parameter 
R distributed around R = 4 for which there is a zero mass state (see Table®). 

A few comments are worth making 

• Referring to the Table, one sees that the average distance Am 8 = 
m i+ i — rrii with i > 2 between two consecutive states decreases very 
slowly along a tower and is of the order, 

Ami w 0.25 TeV . (31) 



• As a result of our numerical computations (k = 1) performed for neigh- 
bouring values of kR, namely 

kR = 6.3 + A(kR) , (32) 

we have found that the corresponding mass eigenvalues in a given tower 
can be deduced precisely by 

m m = m [6.3] e -2.A(S«) at fixed ^ (33) 
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The physical masses of the tower in the s = so brane, i.e. mf^j, are 
obtained from (126]) (taking into account (|27|) ) by 



™£? = ™5 31e ^ ( 34 ) 

and hence, more generally 

m g = m! 6 3l e^"°- 2 - A (^) . (35) 

This is true at k fixed. Remember that the consequences of a change 
of k (at fixed kR) can be obtained from the above rescaling properties 

(EES). 

The differences between the masses in the towers are exponentially 
sensitive to the parameter k(so — 2nA(R)) while it is multiplicatively 
sensitive to k. Should a Kaluza-Klein tower be discovered an approxi- 
mate value of k(so — 2nA(R)) could be deduced. 

A particular attention has to be drawn on the first state. Following the 
path in the R parameter space, going up from a large negative value of 
R, one sees that the first mass eigenvalue mi decreases faster than m% 
and gets equal to zero as one reaches the surface R — 4 (k — 1) of the 
zero mass constraint (see Table©). At that point, m 2 — mi is about 
twice the average Am 8 value. Once the zero mass surface is passed, the 
first mass eigenvalue disappears and a tachyon state develops with h 
increasing rapidly. 

The mass m; increases rather slowly when R decreases toward — oo and 
becomes equal to the mass m i+ i corresponding to R = +oo, exhibiting 
continuity of the masses as functions of R. 



Case Al mass eigenvalues 

For the Case Al, the mass equation (see Table (Tjl])) is invariant under the 
rescaling analogous to fl2"g|) . namely [R — > R/X, k — > Xk, a.\ — > di, a 2 — > 
a 2 /A, a 4 — > a 4 , m — ► Am], and hence 

mi(k, ai, 6i2i 6i±) = — (xk, aci, (36) 
mi(k, di\, a.2, 014) = kmi (l, dt\, ka^-, 0.4) ■ (37) 
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The ten lowest mass eigenvalues in the Kaluza-Klein towers are given 
in Table([7]) for an arbitrary choice of the boundary parameters ai = 0.7, 
a4 = 6.6286 and for a set of values of a 2 including the value d 2 = (ai«4 — 
l)/(4aifc) = 1.3 where the path in parameter space crosses the zero mass 
surface. The value a 2 = is not only excluded but appears as a singular 
point. It is then convenient to vary the values of a 2 from zero to +00 and 
then from —00 back to zero. It should be noted that the passage through 
the value a 2 = ±00 is smooth. 

Apart from a few small differences, the main structure of the mass towers 
is essentially the same as for the Case A4 above. The mass eigenvalues and 
the physical masses corresponding to other values of R and s (for k fixed) 
can again be deduced by the same formulae (I3"3"|) - (j3"5|) . A change of k follows 
from the rescaling equation (|30"1) . 

Moreover, comparing the towers which have a zero mass state at their 
bottom (line R = 4 for Case A4 in Table® and line a 2 = 1.3 for Case Al 
in Table©), we note that all the masses in these two towers are practically 
identical. This holds for any kR in the physically allowed range. 

General comments about Case A3 mass eigenvalues 

Summarizing 

• The rescaling [R —> R/X, k 
leads to the formulae 

m i (k,p 1 ,p 2 ) 
m i (k,p 1 ,p 2 ) 

• The approximate parameter condition for the existence of a zero mass 
is (see Table ®) 

p 1 (p 2 -4k) = Q . (40) 

For each of the two solutions p~\ = or p 2 = 4k, there indeed exists a 
zero mass state. 

• For p 2 = 4 (Table®) the mass tower is identical to the mass towers 
with a zero mass state in the cases Al and A4 (Tables® and ®), and 
this independently of the value of p\. 



Xk, pi — > Xpi, p 2 — > Xp 2 , m — > Am], 



1 

A 

km 



mi (Xk, Xpi, Xp 2 ) 
Pi P2 



1 \ 1 j 1 



k 1 k 



(38) 
(39) 
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• For pi = 0, besides the zero mass the other masses depend on the value 
of pi- When p 2 moves toward 4, the lowest non zero mass in the tower 
converges also to zero. 

• When p2 is fixed to a given value, all the towers corresponding to any 
value of pi are identical, including the tachyon if it exists (i.e. for 
p2 > 4). There is of course an extra mass zero for pi = 0. However, 
for pi close to zero, neither a small mass particle nor a small h tachyon 
appears. 

5 Physical Considerations. Probability den- 
sities 

In the context of a given boundary case, once all the parameters are fixed and 
the mass eigenvalue tower is determined, there exists a unique field <pn\s) 
for each mass eigenvalue leading to a normalized probability density field 
distribution D n (s) along the fifth dimension ([5]) 



It is convenient to parametrize the s range [0, 2irR] by the reduced variable 
x defined by 



with range [0, 1]. 

Case Al field probability densities 

We consider the tower labeled by 6:2 = 1-3 in Table([7]). The logarithm of 
the normalized probability density for the three mass eigenvalues mi = 0, 
m.2 = 0.501 TeV and m 5 = 1.275 TeV are given in Figure([T|), (T5]) and ([3]) 
respectively as functions of x. The general trends are as follows: 

• The probability density is a fast varying function of x. In a large part 
of the domain the logarithm increases or decreases linearly. 

• A general pattern emerges. For rrii with i even, the probability density 
presents around x = 0.5 a very steep dip down to zero as a result of a 
brutal but continuous change of sign of the field at that point. 




(41) 
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• Moreover, the mass m; presents i—1 probability dips in the high x 
region 0.95 < x < 1. 

Relative mass eigenstate probabilities for given xq 

In a brane supposed to be located at a certain fixed x (s = 27iRx ), it is 
directly possible to compare the probabilities of the different mass eigenstates 
in a given tower. Neglecting dynamical and kinematical effects related to the 
production in the available phase space, these probabilities given in Table© 
would account for the rate of appearance of the mass eigenvalue states to an 
observer sitting at this x$. Note however that the physical masses as seen by 
this observer (at Xq 7^ 0) are not the eigenvalue masses but vary with the Sq 
in agreement with (1551) . 

The ratios of probabilities densities are given in Table© for the first 
ten mass eigenstates and selected values of x in its range [0,1], arbitrarily 
normalizing to one the highest probability among the ten first masses con- 
sidered. 

• For x values outside the region where dips in the probability density 
appear i.e. in the two regions < x < 0.40 and 0.55 < x < 0.90, 
the relative probabilities are very weakly dependent on x. In the first 
region < x < 0.40, the probability for m\ dominates whereas in the 
second region 0.55 < x < 0.95, it dominates for mi . 

• In the dip regions 0.45 < x < 0.5 and 0.90 < x < 1, a chaotic behavior 
shows up. Small variations of x may imply large fluctuations of the 
probabilities. 

Relative probabilities in a given physical mass tower as a function 

Of Xq 

In Table (fTDl . we have taken into account the increase of R dictated by (I35I) 
(R = 6.3/ (1 — xq)) as requested to keep the low lying physical masses exactly 
unchanged as one increases Xq starting from xq = where mass eigenvalues 
and physical masses coincide. We have limited ourselves to R < 50. Indeed, 
larger values of R obviously spoil the underlying philosophy that all the 
reduced parameters have to be of order one. The similarity between the 
results presented in the two Tables Qj and ffTUj) should be noted though the 
physical interpretation is widely different. 
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Case A4 field probability densities 

Summarizing 

• Figures (j3J) and (jSJ) show the logarithm of the normalized probability 
density as a function of x for R = 4, respectively for two mass eigen- 
values chosen for illustration mi = and = 0.766 TeV. 

• All the probabilities are seen to increase very fast with x since the 
logarithm is essentially a linear function for 0.1 < x < 0.9. 

• Due to the boundary condition 0(0) = (Table ([T])) there is a sharp dip 
for s — 0. 

• In the high x region (0.9 < x < 1) the probability for the mass rrii 
exhibits i—1 sharp dips corresponding to zeros in the field. 

• In the whole region 0.1 < x < 0.9 the relative probabilities in a tower 
are very close to those of Case Al for x between 0.5 and 0.9 (Table 

• In the two extreme regions (where there are dips) the relative proba- 
bilities exhibit a chaotic behavior as in the dip regions for Al. 

• The same considerations hold as for the case Al regarding the requested 
readjustment of R to keep the physical masses unchanged when the 
position xo is changed. 

6 Conclusions 

Inspired by the warped five-dimensional scenario of Randall and Sundrum 
and restricting to the case of a real massless scalar field supposed to propagate 
in the bulk, we have developed a similar warped model, keeping all the basic 
parameters adjusted in terms of the Planck mass as the only dimensionful 
scale. This, in the end, solves the mass scale hierarchy problem. 

We have concentrated primarily on a careful study of the hermiticity 
(symmetry, self-adjointness) and commutativity of the operators susceptible 
to be used to validly establish the Kaluza-Klein reduction equations in the 
five dimensional warped space. Postulating that the fifth extra dimension s is 
compactified on a strip < s < 2nR and that the metric has no dicontinuity, 
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we have enumerated all the allowed boundary conditions. From them we 
have deduced all the Kaluza-Klein towers mass equations providing the mass 
eigenvalues, as well as the tachyon mass equations. 

We have discussed how these mass eigenstates show up with physical 
masses depending on the location of our brane on the s axis. 

As an illustration, we have carried on some numerical computations for 
the three sets of boundary conditions Al, A4 and A3 in order to visualize the 
structure of the towers and to investigate their main properties. The other 
cases can be studied along the same lines. 

The structure of the eigenvalue towers depends generally in a sensitive 
way on the value of the basic parameter kR and to a smaller extend on the 
boundary parameters and on the reduced parameter k. Apart from small 
differences, the main structure is the same in the three Cases considered. 

One notices that the first tower masses are of the order TeV for kR around 
6.3 as expected for solving the mass hierarchy problem. 

In general, for specific values of the boundary parameters there exists a 
zero mass state. For parameters close to these values the first mass state in 
the tower is often either a particle with a small mass or a tachyon. 

One observes that the mass spacing (discarding the first mass state) is 
very stable within a given tower and is exponentially sensitive to the value 
of k(s — 2nR) where s is the position of our brane. Hence should a Kaluza- 
Klein tower be observed experimentally, a good estimation of this basic pa- 
rameter would result. 

The normalized field probability density for any physical mass state in a 
tower can easily be computed for fixed values of the boundary parameters as 
a function of k, R and xq. Neglecting dynamical and kinematical effects the 
ratios of the probabilities among the first masses in a given tower evaluated 
at Xq would express the relative intensities of the eventually observed mass 
peaks. 
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A Appendix: The Bessel Functions. Nota- 
tions 

The Bessel function J p , as well as Y p , satisfies 

y 2 d 2 y +yd y + (y 2 -p 2 )y p (y) = 

U(v)+p— - J P -M = 

J p (y) - 2(p - i)i=lM + J p _ 2 (y) = 0. (43) 



A useful identity is 



J 2 (y)Y 1 (y) - Ji(y)y 2 (y) - — = . (44) 

vry 



The modified Bessel function J p , as well as K p , satisfies 

(y 2 d 2 y + yd y - (y 2 + p 2 ))^) = (45) 
and (note the sign differences) 

dyiAy)+p Ij y--i P -i(y) = o 

d y K p {y) + + K p ^{y) = 

I p ( y ) + 2(p-l) I -^--I p _ 2 (y) = 
^(y)-2(p-l)^=iM-X p _ 2 (y) = 0. (46) 

y 

The corresponding useful identity is 

I 2 (y)K 1 (y) + I 1 (y)K 2 (y)-^0. (47) 
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B Appendix: Singularities 



In the main part of the paper, we have analyzed the case of a warped space ([I]) 
induced by a cosmological constant §2§ with a fifth dimension compactified 
to a strip < s < 2tcR. We made the choice of a positive k everywhere in 
the space and postulated that there was no singularity 

First, the k negative case is simply related to the positive case by the 
exchange s <-> — s+2irR which hence makes the sign of k an arbitrary choice. 

However, it may be assumed that in some region of s the constant k 
is positive and in another region it is negative. With either k positive or 
negative there must be at least a singular point s s where a transition in the 
metric occurs. Writing 

for s < s s dS 2 = Ce 2k{s ' Sa) dx^ - ds 2 

for s = s s dS 2 = C dx^dx^ — ds 2 

for s > s s dS 2 = Ce-^-^dx^ - ds 2 (48) 

the metric is continuous at s = s s as it should but its first derivative has 
a discontinuity ±4fcC and its second derivative a 5-function behavior. In 
principle, there could be any finite number of such singularities. 

If there is an even number of singularities (2m, m > 1), the strip compact- 
ification can be transformed to an orbifold compactification by identifying the 
edges and 2irR closing the strip to a circle, allowing periodic or antiperiodic 
conditions. Without loss of generality, one may chose these 2m singularities 
to be located at 0, si, s 2 , . . . , s 2?n _i with < si < s 2 < . . . < s 2m -i < 2nR. 
A necessary condition for the closure is 

m m—1 
i=l i=l 

Indeed the total length of the region where k is positive must be equal to the 
total length where k is negative and must thus be equal to one half of the 
length 2nR of the circle. 

We expect to come back to the problem of the warped Kaluza-Klein 
towers with singularities in a forthcoming paper. This change of metric has 
a direct impact on the establishment of the boundary conditions and the 
treatment of the Kaluza-Klein equations. 



19 



References 



[1] Grard, F., Nuyts, J., Phys.Rev. D 74 , 124013 (2006), [hep^th/0607246 



[2] Arkhani-Ahmed, N., Dimopoulos, S., Dvali, G., Phys. Lett., B429, 263 
(1998), |hep^ph/9803315[ SLAC-PUB-7769, SU-ITP-98/13 

Antoniadis, I., Arkani-Hamed, N., Dimopoulos, S., Dvali, G.R., 
Phys. Lett, B436, 257 (1998). [hep^ph/9804398[ SLAC-PUB-7801, SU- 
ITP-98-28, CPTH-S608-0498 

[3] Kaluza, T., Sitzungsber. Preuss. Akad. Wiss. Berlin. (Math. Phys.), 966- 
972 (1921). Klein, O., Z. Phys. 37, 895-906 (1926). 

[4] Asorey, M., Garcia Alvarez, D., Munoz-Castaneda, J.M., 



hep-th/ 0604089 



[5] Randall, L., S undrum, R., Phys. Rev. Lett. 83 , 3370 (1999) 
hep-ph/990522ll Physical Review Letters 83 , 4690 (1999), 
hep-t h/99060"64| 



[6] Bender, CM., Boettcher, S., Phys. Rev. Lett. 80 , 5243 (1998), 



|physics/971 2001 



[7] Fairlie, D.B., Nuyts, J., J. Phys. A38, 3611-3624 (2005),[hep^th/0412148 



20 



Table 1: Two boundary conditions 



Two Boundary Conditions 


Case 


Boundary Conditions 


Reduced Parameters 


Al 


4>{2tvR) = e 47r/crt (ai0(O) + a 2 d s <f>{0)) 
d 3 <fr(2irR) = e^ kR («iHi=i0( O ) + a 4 9 s 0(O)) 


ai = ai , ct2 = a2/Mpi ^ 
az4 = 0C4 


A2 


4>{2ttR) = e 47rfcK ai0(O) 

ds^R) = e^ kR (a 3 0(O) + ^d s 4>(0)) 


ai = ai 7^ 

Q?3 = «3 M Pi 


A3 


&0(O) = pi 0(0) 
ds4>i2irR) = P24>{2tvR) 


pi = pi Mpi 
f>2 = p2 Mpi 


A4 


m = 

ds4>{2iTR) = k4>(2ttR) 


K = K Mpi 


A5 


0(2ttR) = 
<9 S 0(O) = C<H0) 


( = (M Pl 


A6 


0(0) = 
0(2ttR) = 





Table 2: Three boundary conditions 



Three Boundary Conditions 


Case 


Boundary Conditions 


Reduced Parameters 


Bl 


0(2ttR) = e i7rkH \ 1 <j)(0) 
ds4>(0) = A 2 0(O) 
a s 0(27rR) = e 47rfc «A 3 0(O) 


Ai = Ai 

<^2 = A2 Mpi 
A3 = A 3 M Pi 


B2 


0(0) = 

4>(2itR) = e 47rfci Vid s 0(O) 
9 s 0(27ri?) = e 4 **^ 9*0(0) 


pi = pi/Mp; 
P2 = P2 


B3 


0(0) = 

a s 0(o) = o 

a 3 0(27rii) = ucP(2ttR) 


v = vMpi 


B4 


0(0) = 
0(2ttR) = 
d s 0(0) = 
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Table 3: Two boundary conditions : zero mass state constraints (— > : ap- 
proximate relation resulting from e 4k7rR being very large) 



Parameter constraints for zero mass states 


Case 


Parameter constraint 


uo 5 To relation 


Al 


-8 e ilTkR a 2 k 

+ («ia4 + 4a2Ct4fc — 1) = 


g^fcH (e 4 '^ - ai - 4a 2 fc) <t = (e 47TkH ai - l) ro 


— » Aaict'jk — aia^ + 1 S3 




A2 


-8 e 4 *" fcfl aiifc 

+ (aiQ3 + 4fc) = 




— > 4a^fc — o?3 Ri 


-V CT() Si 


A3 


e 8^ fcKpi (4fe _ p2) + p2 (pi _ 4fc) = o 


(4fc - pi) (TO = pi TO 


— ► Pi (p2 - 4fe) S3 


— ► (4fe - pi) (TO = pi T() 


A4 


gSvrfcK ( 4 fc _ K ) + K = 


TO = "TO 


-+ K - 4fc S3 


— ► (T{) = —TO 


A5 


e SvrfcK£ + ( 4fc _ =Q 


e«™ K a = -T 


-» C~o 


->■ (TO S3 


A6 


No zero mass state 


(TO — To — 
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Table 4: Mass tower equations for two boundary conditions 





Notations: E = e 2wkR , F = 


. m TTi „2nkR m 

- k , i-l — e fc 


Case 


Mass equation 


cr, r relation 


Al 


a 2 mE[(J 2 (F )Y 1 (F 2 ) - J 1 (F 2 )Y 2 (F )) «i 
+a 2 m(J 1 (F )Y 1 (F 2 ) - Ji(F 2 )Y 1 (F ))] 

+ [(ai« 4 - I) {J 2 {F 2 )Y 2 {F ) - J 2 (F )Y 2 (F 2 )) 
+a 2 a 4 m(J 2 (F 2 )Y 1 (F ) - J 1 (F )Y 2 (F 2 )) 


(J 2 (F )ai + ,/i(F )a 2 m - J 2 (F 2 )) CT 

+ (Y 2 (F )a 1 + Yi(F )a 2 m - y 2 (F 2 ))r = 


A2 


mEal [J 2 (F )Y 1 (F 2 ) - J 1 (F 2 )Y 2 (F )] 
+ [(./ 2 (F 2 )Y 2 (F ) - J 2 (F )Y 2 (F 2 )) ai a 3 

+m(J 2 (F 2 )Y 1 (F ) - Ji(F )y 2 (F 2 ))] = 


(J 2 (F )ai - ,/ 2 (F 2 )) CT 

+ (Y 2 (F )a 1 - y 2 (F 2 ))r = 






A3 


mF[ P1 (Ji(F 2 )Y 2 (F ) - J 2 (F )Y 1 (F 2 )) 
+m(J 1 (F )Y 1 (F 2 ) - J^F^Y^Fg))] 

+P 2 [P1 (J 2 (F )Y 2 (F 2 ) - J 2 (F 2 )Y 2 (F )) 

+m,(J 2 (F 2 )Y 1 (F ) - Ji(F )y 2 (F 2 ))] = 


(Jl(Fo)m- J 2 (F )pi)a 

+ (Y 1 (F )m - Y 2 (F )pi)t = 


A4 


mi5(Ji(F 2 )Y 2 (F ) - ./atFoJ^fFa)) 

+ K (j 2 (F )y 2 (F 2 ) - J 2 (F 2 )y 2 (F )) = o 


J2(F )<t + Y 2 (F )t = 


A5 


C(.7 2 (F 2 )y 2 (F ) - J 2 {F )Y 2 (F 2 )) 

+m(J 1 {F )Y 2 (F 2 ) - J 2 (F 2 )y 1 (F )) = 


J2(F 2 )<r + Y 2 (F 2 )t = 


A6 


./2(F 2 )r 2 (F ) - J 2 {Fo)Y 2 {F 2 ) = 


J 2 (F )a + Y 2 (F )t = 



23 



Table 5: Tachyon equations for two boundary conditions 



Notations: E = e 27lkR , F = £ , F 2 = e 27rkn \ (m 2 = -h 2 ) 


Case 


Mass equation 


a, r relation 


Al 


a 2 hE[(-I 2 (F )K 1 (F 2 ) - h(F 2 )K 2 (F )) ai 
+a 2 h(-I 1 (F )K 1 (F 2 ) + /l(F 2 )Ki(F ))J 

+ [(ai« 4 - l)(l2(F 2 )K 2 (F ) - I 2 (F )K 2 (F 2 )) 
+a 2 a 4 h(-I 2 (F 2 )K 1 (F ) - I 1 (F )K 2 (F 2 )) 
+ 2a 2 k] = 


(I 2 (F ) ai + Ii(F a )a 2 h - I 2 (F 2 ))<j 

+ (K 2 (F )a 1 - K 1 (F )a 2 h - K 2 (F 2 )) t = 


A2 


a 2 hE[(-I 2 (F )K 1 (F 2 ) - h(F 2 )K 2 (F )) 
+ [(I 2 {F 2 )K 2 (F ) - I 2 (F )K 2 (F 2 )) aia 3 

+h (-72(^2)^1(^0) - Il(F )K 2 (F 2 )) 

+2 ai k] = 


(I 2 (F ) ai - I 2 {F 2 ))a 

+ (K 2 (F ) ai - K 2 (F 2 ))t = 


A3 


hE[ P1 (h(F 2 )K 2 (F ) + I 2 (F )K 1 (F 2 )) 

+h(-I 1 (F )K 1 (F 2 ) + I 1 (F 2 )K 1 (F ))] 

+P2[P1 (72(^0)^2(^2) - /2(F 2 )*r 2 (Fo)) 

-h(I 2 (F 2 )K 1 (F ) + I 1 {F )K 2 (F 2 ))} = 


(I 1 (F )h - I 2 (F ) pl )a 

- (K 1 (F )h + K 2 (F ) P1 ) t = 


A4 


hE(I 1 (F 2 )K 2 (F ) + I 2 (F )K 1 (F 2 )) 

+k(I 2 (F )K 2 (F 2 ) - I 2 (F 2 )K 2 (F )) = 


7 2 (F )<t + K 2 (F )t = 


A5 


t(l 2 (F 2 )K 2 (F ) - I 2 (F )K 2 (F 2 )) 

+h(l 1 (F )K 2 (F 2 ) + I 2 (F 2 )K 1 (F a )) = 


7 2 (F 2 )<t + K 2 (F 2 )t = 


A6 


I 2 (F 2 )K 2 (F ) - I 2 (F )K 2 (F 2 ) = 


7 2 (F )<t + K 2 (F„)t = 
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Table 6: 



Towers of mass eigenvalues for the Case A4, masses are in TeV, kR = 6.3, k = 1 



K 


h 


mi 


ni2 


m 3 


ir»4 


m 5 


m 6 


mi 


m 8 


m 9 


mio 


- 100.0 




0.4 


0.655 


0.9048 


1.152 


1.398 


1.644 


1.889 


2.134 


2.379 


2.624 


- 12.0 




0.3769 


0.62 


0.8599 


1.1 


1.34 


1.582 


1.825 


2.068 


2.312 


2.557 


- 8.0 




0.3672 


0.6072 


0.846 


1.086 


1.328 


1.571 


1.814 


2.058 


2.303 


2.548 


- 4.0 




0.349 


0.5866 


0.8265 


1.069 


1.312 


1.557 


1.802 


2.047 


2.293 


2.539 


0.0 




0.301 


0.5512 


0.7993 


1.047 


1.294 


1.542 


1.789 


2.036 


2.283 


2.529 


3.9 




0.0605 


0.5025 


0.768 


1.024 


1.276 


1.526 


1.775 


2.024 


2.272 


2.52 


3.99 




0.0192 


0.5014 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.272 


2.52 


3.999 




0.00609 


0.5016 


0.766 


1.023 


1.275 


1.525 


1.775 


2.023 


2.272 


2.52 


4.0 







0.5013 


0.766 


1.023 


1.275 


1.525 


1.775 


2.023 


2.272 


2.52 


4.001 


0.00609 




0.5016 


0.766 


1.023 


1.275 


1.525 


1.775 


2.023 


2.272 


2.52 


4.01 


0.01926 




0.5016 


0.766 


1.023 


1.275 


1.525 


1.775 


2.023 


2.272 


2.52 


4.1 


0.0612 




0.501 


0.766 


1.022 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


8.0 


0.4826 




0.4624 


0.7380 


1.0 


1.256 


1.51 


1.761 


2.011 


2.261 


2.51 


12.0 


0.8091 




0.4428 


0.718 


0.9819 


1.241 


1.496 


1.749 


2.0 


2.251 


2.501 


16.0 


1.128 




0.4327 


0.7058 


0.9689 


1.228 


1.484 


1.737 


1.99 


2.241 


2.492 


20.0 


1.445 




0.426 


0.6969 


0.959 


1.218 


1.474 


1.728 


1.981 


2.233 


2.484 


100.0 


7.737 




0.408 


0.6689 


0.9227 


1.175 


1.426 


1.676 


1.926 


2.176 


2.426 



Table 7: 



Towers of mass eigenvalues for the Case Al, masses are in TeV, kR = 6.3, k = 1, a\ = 0.7, «4 = 6.6286 



Q?2 


h 


mi 


rrt 2 


m 3 


r«4 


m 5 


m 6 


mi 


TOg 


mg 


mm 


0. 


oo 




0.4035 


0.6613 


0.9129 


1.162 


1.411 


1.659 


1.907 


2.154 


2.402 


0.02 


20.31 




0.4051 


0.6639 


0.9165 


1.167 


1.416 


1.665 


1.914 


2.163 


2.411 


0.2 


1.918 




0.4204 


0.6884 


0.949 


1.206 


1.462 


1.716 


1.969 


2.222 


2.473 


1. 


0.2269 




0.4873 


0.7575 


1.015 


1.269 


1.521 


1.771 


2.019 


2.268 


2.516 


1.29 


0.0340 




0.5012 


0.766 


1.022 


1.274 


1.525 


1.774 


2.023 


2.271 


2.519 


1.299 


0.0107 




0.5013 


0.7669 


1.022 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


1.3 







0.5013 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


1.301 




0.01064 


0.5013 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


1.31 




0.03359 


0.5017 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


1.6 




0.159 


0.5112 


0.773 


1.027 


1.278 


1.528 


1.777 


2.025 


2.273 


2.521 


2.6 




0.24 


0.5269 


0.7833 


1.035 


1.284 


1.533 


1.781 


2.029 


2.277 


2.524 


40.0 




0.2984 


0.5497 


0.7983 


1.046 


1.293 


1.541 


1.788 


2.035 


2.282 


2.529 


100.0 




0.3 


0.5511 


0.7989 


1.046 


1.294 


1.541 


1.788 


2.035 


2.282 


2.529 


±oo 




0.3010 


0.5512 


0.7993 


1.047 


1.294 


1.541 


1.788 


2.035 


2.282 


2.529 


- 100.0 




0.3021 


0.5518 


0.7997 


1.047 


1.294 


1.541 


1.788 


2.035 


2.282 


2.529 


- 40.0 




0.3037 


0.5527 


0.8011 


1.048 


1.295 


1.542 


1.789 


2.035 


2.282 


2.529 


- 2.0 




0.3375 


0.5763 


0.8178 


1.061 


1306. 


1.551 


1.797 


2.043 


2.289 


2.535 


- 0.001 




0.4034 


0.6612 


0.9127 


1.162 


1.411 


1.659 


1.906 


2.154 


2.402 


2.648 
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Table 8: 



Towers of mass eigenvalues for the Case A3, masses are in TeV, kR = 6.3, k = 1 



P2 


h 


mi 


rti2 


7713 


IT14 


ms 


me 


mj 


mg 




TO10 


- 4000.0 




0.4034 


0.6611 


0.9127 


1.162 


1.411 


1.659 


1.906 


2.154 


2.401 


2.648 


- 400.0 




0.4025 


0.6596 


0.9106 


1.160 


1.407 


1.655 


1.902 


2.149 


2.396 


2.642 


- 40.0 




0.3940 


0.6459 


0.8919 


1.136 


1.380 


1.623 


1.866 


2.109 


2.352 


2.596 


- 4.0 




0.3487 


0.5866 


0.8265 


1.069 


1.312 


1.557 


1.802 


2.047 


2.293 


2.539 


0.0 




0.3010 


0.5512 


0.7993 


1.047 


1.294 


1.541 


1.788 


2.035 


2.282 


2.529 


3.0 




0.1807 


0.5139 


0.7750 


1.029 


1.279 


1.529 


1.778 


2.026 


2.274 


2.522 


3.9 




0.0605 


0.5025 


0.766 


1.023 


1.275 


1.525 


1.775 


2.023 


2.271 


2.519 


3.99 




0.0192 


0.5014 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


4.0 




0. 


0.5012 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


4.01 


0.0193 




0.5011 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


4.1 


0.0612 




0.5000 


0.766 


1.022 


1.274 


1.524 


1.774 


2.023 


2.271 


2.519 


5.0 


0.2046 




0.4895 


0.7590 


1.017 


1.270 


1.521 


1.771 


2.020 


2.268 


2.517 


8.0 


0.4826 




0.4624 


0.7380 


1.0 


1.256 


1.509 


1.761 


2.011 


2.260 


2.509 


12.0 


0.8092 




0.4425 


0.7180 


0.9818 


1.240 


1.495 


1.748 


2.0 


2.250 


2.500 


40.0 


3.021 




0.4142 


0.6786 


0.9365 


1.192 


1.446 


1.699 


1.952 


2.204 


2.455 


400.0 


31.31 




0.4045 


0.6629 


0.9152 


1.165 


1.415 


1.663 


1.912 


2.160 


2.408 


4000.0 


314.1 




0.4036 


0.6614 


0.9131 


1.163 


1.411 


1.659 


1.907 


2.155 


2.402 
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Table 9: The relative probability densities of the ten first mass eigenvalues 
rrii o (in TeV) for given x = s/ (2ttR) in the tower corresponding to the Case 
Al with kR = 6.3, k = 1, 6t\ = 0.7, «2 = 1.3, = 6.6286. The highest 
probability among these ten masses is normalized to exactly 1 and labeled 
as such. Note that the physical masses depend on s , hence on x, and are 
given in terms of the mass eigenvalues by m iiS0 = rriifie 27TxkR in agreement 
with Eq.flMD. 



Mass 
eigenvalues 


W»1,0 


"*2,0 


m 3 ,o 


"M,o 


"15,0 


"»6,0 


"*7,0 


m8,o 


"»9,0 


"1-10,0 





0.5013 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


X 


Relative Probabilities 





1 


0.27 


0.19 


0.14 


0.11 


0.096 


0.083 


0.073 


0.065 


0.059 


0.4 


1 


0.27 


0.19 


0.14 


0.11 


0.096 


0.083 


0.073 


0.065 


0.059 


0.45 


1 


0.27 


0.19 


0.14 


0.12 


0.084 


0.10 


0.055 


0.090 


0.034 


0.46 


1 


0.26 


0.21 


0.11 


0.17 


0.044 


0.18 


0.0085 


0.23 


0.0011 


0.466 


1 


0.24 


0.25 


0.067 


0.27 


0.0024 


0.41 


0.035 


0.69 


0.22 


0.467 


1 


0.24 


0.26 


0.057 


0.30 


0.000021 


0.49 


0.071 


0.86 


0.35 


0.4676 


1 


0.24 


0.27 


0.051 


0.33 


0.00066 


0.55 


0.10 


0.99 


0.45 


0.4676334 


1 


0.24 


0.27 


0.050 


0.33 


0.00075 


0.55 


0.10 


1.0 


0.46 


0.4676335 


1.0 


0.24 


0.27 


0.050 


0.33 


0.00075 


0.55 


0.10 


1 


0.46 


0.4677 


0.98 


0.23 


0.26 


0.049 


0.33 


0.00095 


0.55 


0.11 


1 


0.46 


0.468 


0.92 


0.21 


0.25 


0.042 


0.32 


0.0021 


0.54 


0.12 


1 


0.49 


0.47 


0.56 


0.12 


0.17 


0.013 


0.26 


0.018 


0.50 


0.19 


1 


0.67 


0.475 


0.14 


0.022 


0.068 


0.0017 


0.16 


0.081 


0.42 


0.36 


0.96 


1 


0.48 


0.026 


0.0017 


0.027 


0.012 


0.10 


0.11 


0.31 


0.39 


0.78 


1 


0.5 


0.00012 


0.0013 


0.0092 


0.027 


0.067 


0.13 


0.25 


0.41 


0.66 


1 


0.9 


0.000022 


0.0017 


0.0092 


0.029 


0.069 


0.14 


0.25 


0.42 


0.67 


1 


0.92 


0.000028 


0.0022 


0.011 


0.035 


0.083 


0.16 


0.29 


0.46 


0.70 


1 


0.94 


0.00011 


0.0081 


0.040 


0.11 


0.23 


0.40 


0.59 


0.78 


0.93 


1 


0.95 


0.00052 


0.036 


0.16 


0.40 


0.69 


0.93 


1 


0.86 


0.56 


0.24 


0.97 


0.012 


0.51 


1 


0.58 


0.017 


0.21 


0.40 


0.095 


0.044 


0.26 


0.99 


0.71 


0.85 


1 


0.084 


0.24 


0.49 


0.065 


0.13 


0.33 


0.060 


0.995 


1 


0.0075 


0.081 


0.25 


0.31 


0.21 


0.064 


0.000017 


0.048 


0.13 


0.999 


1 


0.26 


0.16 


0.11 


0.073 


0.049 


0.031 


0.019 


0.010 


0.0044 


1 


1 


0.27 


0.19 


0.14 


0.11 


0.096 


0.083 


0.073 


0.065 


0.059 
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Table 10: The relative probability densities of the ten first physical masses vtii 
(in TeV) for given xq = so/(2tcR) in the tower corresponding to the Case Al, 
k — 1, «i = 0.7, a.2 = 1.3, 0:4 = 6.6286. The highest probability among these 
ten masses is normalized to 1 and labeled as such. When xq is increased, 
R is adjusted according to (1511) . keeping k — 1, so as to retain the physical 
masses unchanged. 



Physical 
masses 


m± 


m 2 


m 3 


1J14 


m 5 


m 6 


my 




m 9 


mio 





0.5013 


0.766 


1.023 


1.275 


1.525 


1.774 


2.023 


2.271 


2.519 


XQ 


if. 


Relative Probabilities 





6.3 


1 


0.28 


0.19 


0.14 


0.11 


0.096 


0.083 


0.073 


0.065 


0.059 


0.44 


11.3 


1 


0.27 


0.18 


0.14 


0.11 


0.095 


0.082 


0.072 


0.064 


0.058 


0.475 


12.0 


1 


0.27 


0.19 


0.13 


0.13 


0.077 


0.11 


0.046 


0.1 


0.025 


0.4793 


12.1 


1 


0.26 


0.21 


0.11 


0.17 


0.041 


0.18 


0.0068 


0.24 


0.0022 


0.4815 


12.15 


1 


0.25 


0.23 


0.081 


0.22 


0.012 


0.3 


0.0062 


0.47 


0.088 


0.4828 


12.18 


1 


0.24 


0.25 


0.059 


0.29 


0.00041 


0.45 


0.056 


0.79 


0.29 


0.4832 


12.19 


1 


0.23 


0.26 


0.051 


0.32 


0.00033 


0.52 


0.093 


0.94 


0.42 


0.4834 


12.195 


0.97 


0.22 


0.26 


0.045 


0.32 


0.0015 


0.55 


0.11 


1 


0.48 


0.4836 


12.2 


0.88 


0.2 


0.24 


0.037 


0.31 


0.0033 


0.54 


0.13 


1 


0.51 


0.4857 


12.25 


0.32 


0.062 


0.11 


0.0012 


0.21 


0.046 


0.47 


0.28 


1 


0.85 


0.49 


12.3 


0.095 


0.013 


0.053 


0.0041 


0.14 


0.09 


0.38 


0.37 


0.9 


1 


0.61 


16.3 


0.000021 


0.0017 


0.0088 


0.028 


0.066 


0.14 


0.25 


0.42 


0.66 


1 


0.76 


26.3 


0.000022 


0.0017 


0.0089 


0.028 


0.067 


0.14 


0.25 


0.42 


0.69 


1 


0.89 


56.3 


0.000021 


0.0016 


0.0081 


0.025 


0.07 


0.12 


0.25 


0.4 


0.6 


1 
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Figure 1: Case Al : The logarithm of the field probability density as a 
function of x = s/ (2nR) for k = 1, kR = 6.3, a\ = 0.7, a 2 = 1.3, a 4 = 6.6286 
and rrii — 



Case A1, kR = 6.3, m<i = 

20 t 




-80 A — — i — — i — — i — — i — — I 
0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 2: Case Al : The logarithm of the field probability density as a 
function of x = s/ (2nR) for k = 1, kR = 6.3, 6l\ = 0.7, a 2 = 1.3, a 4 = 6.6286 
and m 2 = 0.501TeV 
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Figure 3: Case Al : The logarithm of the field probability density as a 
function of x — s/ {2nR) for k — 1, kR = 6.3, ai = 0.7, a 2 = 1-3, a 4 = 6.6286 
and m.5 — 1.275TeV 



Case A1 , kR = 6.3, m = 1 .275 TeV 
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Figure 4: Case A4 : The logarithm of the field probability density as a 
function of x — s/ (2nR) for k — 1, kR = 6.3, R = 4 and m x = 
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Figure 5: Case A4 : The logarithm of the filed probability density as a 
function of x — sj (2nR) for k — 1, kR = 6.3, R = 4 and m 3 = 0.766 TeV 
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